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'. Abstract . In this paper we discuss recent progress in using the Camassa-Holm equations 
to model turbulent flows. The Camassa-Holm equations, given their special geometric 

■ and physical properties, appear particularly well suited for studying turbulent flows. We 
identify the steady solution of the Camassa-Holm equation with the mean flow of the 

! Reynolds equation and compare the results with empirical data for turbulent flows in 



channels and pipes. The data suggests that the constant a version of the Camassa-Holm 



m 
m 

O ■ equations, derived under the assumptions that the fluctuation statistics are isotropic and 

O I homogeneous, holds to order a distance from the boundaries. Near a boundary, these 

Q\ . assumptions are no longer valid and the length scale a is seen to depend on the distance 

, . to the nearest wall. Thus, a turbulent flow is divided into two regions: the constant a 

■ region away from boundaries, and the near wall region. In the near wall region, Reynolds 

O [ number scaling conditions imply that a decreases as Reynolds number increases. Away 

^ , from boundaries, these scaling conditions imply a is independent of Reynolds number. 

^ Given the agreement with empirical and numerical data, our current work indicates that 

^ the Camassa-Holm equations provide a promising theoretical framework from which to 

H ' understand some turbulent flows. 
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Introduction . Laminar Poiseuille flow occurs when a fluid in a straight channel, or pipe, is 
driven by a constant upstream pressure gradient, yielding a symmetric parabolic stream- 
wise velocity proflle. In turbulent states, the mean streamwise velocity remains symmetric, 
but is flattened in the center because of the increase of the velocity fluctuation. Although a 
lot of research has been carried out for turbulent channel flow,^~^ accurate measurement of 
the mean velocity and the Reynolds stress profiles, in particular for flows at high Reynolds 
numbers, is still an experimental challenge. However, in the case of pipe flow, recent ex- 
periments for measuring the mean velocity proflle have been successfully performed for 
moderate to high Reynolds numbers by Zagarola^ The fundamental understanding of how 
these proflles change as functions of Reynolds number, however, seems to be still missing. 

In wall bounded flows it is customary to deflne a characteristic velocity tt* and wall- 
stress Reynolds number Rq by tt* = y^\To\/p and Rq = du^/v, where tq is the boundary 
shear stress. We take the density p to be unity, v is the molecular viscosity of the fluid, and 
d is a characteristic macrolength. For instance, for channel flow d is the channel half-width, 
and for pipe flow d is the pipe radius. Based on experimental observation and numerical 
simulation, a piecewise expression of the mean velocity across the channel or the pipe 
has been commonly accepted,* for which the nondimensional mean streamwise velocity, 
= t//tt*, is assumed to depend onr] = u^z/v and have three types of behavior depending 
on the distance away from the wall boundary, z : a viscous sublayer, in which (f) ^ rj; the 
von Karman-Prandtl logarithmic "law of the wall," in which (f){rj) = K~^ln?7 -|- A, where 
K ~ 0.41 and A ~ 5.5; and a power law region, in which ~ 77^, < p < 1. Alternatively, 
a single curve fltting over the whole region may be proposed (see Panton^). Yet another 
possibility is a family of power laws that flts the data away from the viscous sublayer, and 
has the log law as an envelope, as proposed by Barenblatt et al.^^ 

In this paper (a summary of which was given earlier ■'^■'^), we propose the viscous 
Camassa-Holm equations (VCHE) in (2.14) as a closure approximation for the Reynolds 
equations. The analytic form of our proflles based on the steady VCHE away from the vis- 
cous sublayer, but covering at least 95% of the channel, depends on two free parameters: 
the flux Reynolds number R = dUs^v/t^ (where Ug^y is the streamwise velocity, averaged 
across the channel), and the wall-stress Reynolds number Rq. Due to measurement limi- 
tations most experimental data are contained in this region. Let us remark that we can 
further reduce the parameter dependence to one free parameter by using a drag law for the 
wall friction D ~ Rq/R'^. For the remaining part of the channel, we are unable to solve 
explicitly for the mean proflle without further assumptions, but we do show compatibility 
of the steady VCHE with empirical and numerical velocity proflles in this subregion. The 
VCHE proflles agree well with data obtained from measurements and simulations of turbu- 
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lent channel and pipe flow. For another global approach to turbulent flows in channels and 
pipes displaying good agreement of theoretical mean velocity profiles with experimental 
data see Markus and Smith. 



1. The Euler-Poincare equations and the Euler equations . Consider the Lagrangian com- 
prised of fluid kinetic energy and the volume preservation constraint 
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-Xit,a) 



+ q{X{t, a),t)(detx;(t,a)-l) 



(1.1) 



In (1.1), X(t, a) is the Lagrangian trajectory of the fluid parcel starting at position a at 
time t — 0. Other notation is 



and D{x, t) = (det X'^{t, a))~^ at x = X{t, a) . 



(1.2) 



Moreover, the Jacobian D satisfles the equation 

d 



dt 



D + V ■ (Du) = 



(1.3) 



The extremality conditions for u, where q is viewed as a Lagrange multiplier, are given by 
the Euler-Poincare equation [13] 



(above and throughout we use Einstein's notation for summations) and 

5C 



Sq 







Since 
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the relations (1.3), (1.4), (1.5) yield the Euler equations 

d 



dt 



+ U-V ]u = -Vg, V • « = 



(1.4) 



(1.5) 



Remark. The Euler-Poincare equation (1.4) is equivalent in the Eulerian picture to the 
corresponding Euler-Lagrange equation for fluid parcel trajectories for Lagrangians such 
as (1.1) that are invariant under the right-action of the diffeomorphism group, see Holm 
et al.-*^^"-*^^ and references therein. In what follows, we shall introduce random fluctuations 
into the description of the fluid parcel trajectories in the Lagrangian L in (1.1), take its 
statistical average and use the Euler-Poincare equation (1.1) to derive Eulerian closure 
equations for the corresponding averaged fluid motions. 



2. Averaged Lagrangians and the Camassa-Holm equations . In the presence of random 
fluctuations the Lagrangian trajectory given by X{t, a) has to be augmented with fluctu- 
ations as 

X'^it, a) = X{t, a) + a{X{t, a), t) . (2.1) 

Here a = a{x,t) = a{x,t;Lo) is a random vector field. Thus the Lagrangian L — Liui) 
becomes a random variable 

' + q%X'^{t,a),mei{X^)',{t,a) - . (2.2) 



.(a;) = j daj^ 



In (2.2), we introduce the Eulerian velocity field, 

u'^iy, t) = ^X-(t, a) for y = X-(t, a) , (2.3) 

with X°'(t, a) given in equation (2.1). This is similar to the classical Reynolds decomposi- 
tion of fiuid velocity into its mean and fiuctuating parts. However, this decomposition is 
applied on Lagrangian fiuid parcels, rather than at fixed Eulerian spatial positions. 

Introducing the decomposition (2.3) into the Lagrangian L in (2.2) and changing the 
variables a to x = X{t^ a) yields 

L{uj) = j (ix|y|w"(x + a(x,t),t)|2 + g-(x + a(x,t),t)[det((X<^);oX-i) , 

where D as before is given by (1.2) and satisfies (1.3). Noting that the composition of 
maps X"^ and X gives {X"' o X~^){x, t) = x + a{x, t) we conclude with 

L{uj) = J dxi^^\u''{x + a{x,t),t)\'^ + q''{x + a{x,t),t)[det{I + a'^)-D]^ . (2.4) 

At this stage we make the crucial assumption that a is sufficiently small that the following 
Taylor expansions may be truncated at linear order: 

u^{x + a{x, t), t) ~ u{x, t) + {(t{x, t) ■ W)u{x, t) 

q^ix + a{x, t), t) ~ q{x, t) + {a{x, t) ■ V)q{x, t) (2.5) 
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where 

u{x, t) = {u'^{x + a{x, t),t)) , 

q{x,t) = {q''{x + a{x,t),t)), (2.6) 

and ( ■ ) denotes averaging with respect to the random event co. Thus at this level of 
approximation (2.4) becomes 



1 1 

-\u{x,t)'^\ +u{x,t) ■ {a{x,t) ■ \/u{x,t)) + -\{a{x,t) ■ W)u{x,t)\'^ 



L{uj) = J dxi^D 

+ [q{x, t) + {a{x, t) ■ V)g(x, t)] [det(/ + a'^) - D{x, t)] } . (2.7) 
Therefore the averaged Lagrangian (L) is found to be 

da;|y[|Mp + 2u ■ {{a) ■ V)u+ {aiaj)diU ■ dju]^- 
det(/ + 4)) -D]- D{{a) ■ V)g + {{a det(/ + a',)) ■ V)g} , (2.8) 

where we use the notation di = i = 1, 2, 3. Then the variational derivatives of (L) are 
given by 

= (1 - • iD{a)))u - ^d,{D{a,a,)d,u) 



- (1 + (a) ■ V)q + -[\u\' + 2u ■ ((a) ■ V)u + {a^aj){dju) ■ {d^u)] = -Q (2.9) 
^^^^ = ( det(/ + a',)) -D + V- {{a)D - (adet(/ + a',))) . 



5q 

By stationarity of (L) under variations in q, the last equation in the set (2.9) becomes 

D = (det(/ + a'J) + V • {{a)D) - V • (adet(l + a',)) . 

In order for the mean flow u to be incompressible, one takes D — 1. This imposes the 
condition 

1 = ( det(/ + a'^)) + V • (ct) - V • (cT det(/ + a'J) (2.10) 

on the statistics of the fluctuations. Under this condition, the Euler-Poincare equation 
(1.4) and the equation (1.3) (for (L) instead of L) can be written as 

—v+{u-V)v + VjVuj+VQ = 0, with V • w = , (2.11) 

l/6 



where we define 



V = 



1 d{L) 



(1 - V • {cr))u - di{{aiaj)dju) . 



(2.12) 



D Su 



J D=l 



These equations are shght generahzations of the n-dimensional Camassa-Holm equations. 
The latter correspond to the case where the isotropy conditions 



hold. If moreover the statistics of a are homogeneous, then is constant. Under this 
form the equations (2.11), (2.12) were originally derived. ^^'^^ That derivation generalizes a 
one-dimensional integrable dispersive shallow water model studied in Camassa and Holm^^ 
to n-dimensions and provides the interpretation of a as the typical mean amplitude of the 
fluctuations as in (2.13). 

Remark. The ideal Camassa-Holm equations, or Euler alpha-model, in (2.11) is formally 
the equation for geodesic motion on the diffeomorphism group with respect to the metric 
given by the mean kinetic energy Lagrangian (L) in equation (2.8), which is right-invariant 
under the action of the diffeomorphism group. See Holm et al.^^ for detailed discussions, 
applications and references to Euler-Poincare equations of this type for ideal fluids and 
plasmas. After the original derivation of equation (2.11) in Euclidean space, ^^'^^ Holm 
et al.^^ and ShkoUer^^ generalized it to Riemannian manifolds, discussed its existence 
and uniqueness on a flnite time interval, and amplifled the relation found earlier^^ of this 
equation to the theory of second grade fluids. Additional properties of the Euler equations, 
such as smoothness of the geodesic spray (the Ebin-Marsden theorem) are also known for 
the Euler-CK equations and the limit of zero viscosity for the corresponding viscous Navier- 
Stokes-a equations is known to be a regular limit, even in the presence of boundaries for 
homogeneous (Dirichlet) boundary conditions. ^^'^^ Some of the most interesting solutions 
of the Euler alpha-model could actually leave the diffeomorphism group due to a loss of 
regularity. (This is seen in the one-dimensional Camassa-Holm equation. ^^) Such solutions 
may be interpreted in the sense of generalized flows, as done by Brenier^^ and Shnirelman.^^ 
A functional-analytic study of the Euler alpha-model is made in Marsden et al.^^ 

Adding viscosity. We note that v in (2.12) represents a momentum. Therefore we 
propose that the viscous variant of (2.11) should take the following form, in which the 
viscosity acts to diffuse this momentum. 



(cr) = , {aiaj) = a^Sij 



(2.13) 



d_ 
dt 



v+ {u- V)v + VjVuj = i^Av - VQ , V • w = . 



(2.14) 
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Again, v is given by (2.12). Throughout we will refer to equation (2.14) with definition 
(2.12) as the viscous Camassa-Holm equations (VCHE), or Navier-Stokes alpha-model 
(NS-ck). The standard Navier-Stokes equations are recovered when a is set to zero. The 
VCHE (2.14) in three dimensions possesses global existence and uniqueness, as well as a 
global attractor whose bounds on fractal dimension show cubic scaling with domain size, 
as expected in the Landau theory of three-dimensional turbulence. The proofs of these 
properties of the VCHE, or NS-a model, are given in Foias et al.-'^^-'' 

Since in (2.14), a appears at power up to 2 and we assume |cr| to be small (at least in 
average), the constraint (2.10) can be given a simpler form by using the approximation 

( det(/ + a'^)) - 1 ~ V ■ (a) + {diai ■ 82(72 - d2ai ■ dia2) + 
+ (920-2 • dsas - 93(72 • 920-3) + (930-3 • 9io-i - 9io-3 • 930-1) . 

Then (2.10) becomes (by neglecting the terms of degree > 3 in o-) 

V • ((V • a)a) - V • (o-) ~ (9io-i • ^20-2 - 920-i • 9io-2)+ 

+ (920-2 ■ 930-3 - 930-2 ■ 920-3) + (930-3 ■ 9io-i - 9io-3 ■ 930-1) . (2.15) 

See Gjaja and Holm^^ for the corresponding derivation of equations in the form (2.11) in 
Generalized Lagrangian mean (GLM) theory with (a) = and no viscosity. We note that 
GLM theory provides no closure. 

3. Comicctioii with Coiitiiiumii Mechanics . A mechanical interpretation of these equations 
may be obtained by rewriting the VCHE (2.14) (in the case where {a) = 0, = constant) 
in the equivalent 'constitutive' form 

^ = divT, T= -pl+ 2i^{l - a'^A)D + 2a^D , (3.1) 

with W ■ u = 0,D = (l/2)(Vtt + Vtt^),Q = (l/2)(Vtt - Vtt^), and co-rotational (Jau- 
mann) derivative given hy D = dD/dt + Dft — flD, with d/dt = d/dt + u - V. In this 
form, one recognizes the constitutive relation for VCHE as a variant of the rate-dependent 
incompressible homogeneous fluid of second grade, ^^'^^ whose viscous dissipation, how- 
ever, is modified by the Helmholtz operator (1 — a^A). Thus, the VCHE, or NS-a closure 
model is not only Galilean invariant; it also satisfies the continuum mechanics principles 
of objectivity and material frame indifference. There is a tradition at least since Rivlin^^ 
of using these continuum mechanics principles in modeling turbulence (see also Chorin^^). 
For example, this sort of approach is taken in deriving Reynolds stress algebraic equation 
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models. ' Rate-dependent closure models of mean turbulence have also been obtained by 
the two-scale DIA approach^^ and by the renormalization group methods. 

4. Closure Ansatz . Since VCHE describe mean quantities, we propose to use (2.14) as 
a turbulence closure model and test this ansatz by applying it to turbulent channel and 
pipe flows. For this purpose we also assume that as long as the boundary effects can be 
neglected, the isotropy conditions (2.13) hold. It is also appropriate to recall that the 
Reynolds equations are the averaged Navier-Stokes equations®'^® 

d 

—u + iu ■ V)u = uAu -Vp- (u-u) ■ V(it - u) , V-M = 0, (4.1) 
ot 

where the upper bar denotes the ensemble average, u is the mean flow, p the mean pressure 
and —{{u — u) ■ 'V){u — u) is the divergence of the Reynolds stresses. Our ansatz asserts 
that: 

a) u is approximatively the solution u of the VCHE with the same symmetry and bound- 
ary conditions as u. 

b) The Reynolds stress divergences are given by appropriate terms in the VCHE found 
by matching equations (2.14) and (4.1). 

5. The Reynolds equations for channel flows . For turbulent channel flow (see, e.g., 
Townsend^^), the mean velocity in (4.1) is of the form u — (U (^), 0, O)*'', withp = P{x, y, z) 
and the Reynolds equations (4.1) reduce to 

-vU" + d,{wu) = -da^P , 

d4wz) = -dyP, d,{w^) = -d,P, (5.1) 

where {u^v^wY^ = u — u is the fluctuation of the velocity in the inflnite channel 
{{x,y,z) G R, — d < z < d}. The (1,3) component of the averaged stress tensor 
T — —pi — u<S)U + z/(V?i + (V-u)*^) is given by (T13) = pU'^z) — (wu). At the boundary, 
the velocity components all vanish and one has the stress condition 

Tro = (Ti3)|_^, = ^C/'(.)|_^,, (5.2) 

upon using (wu) = at 2; = ±d. Hence, the Reynolds equations imply {wv){z) = and 
P = Pq — TQx/d — {w'^){z), with integration constant Pq. 

6. The VCHE for channel flows . Passing to the VCHE in the channel, we denote the 
velocity u in (2.14) by U and seek its steady state solutions in the form U= ([/(z), 0, 0)*^ 
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subject to the boundary condition U{±d) = and the symmetry condition U{z) = U{—z). 
In this particular case, the steady VCHE reduces to, 



- p')uy + v{a^U')"' = -da^TT , = -dyTT , = -d^TT , (6.1) 

where = ((t|), (3 = (cxs) and n = tt + J {U - P'U - {a^U')')Udz. 

In accord with the statistical assumptions in the Reynolds equation, we also take the 
statistics of a to be invariant under horizontal translations. As already mentioned above, 
we will suppose that away from the wall, i.e. for \z\ < do with < do < d we have 

a{z) = ao , I3{z) = , (6.2) 



with constants do and ao to be determined later. The following heuristic argument may 
provide some help in understanding this length-scale ao- Clearly a and /3 must depend 
on d, To, z, the only physical quantities present. Dimensional analysis then implies (with 
two suitable functions / and g) that 

" /f«o,^^l. ^=.fflo,^l, (6.3) 



d \ ii, ) d 
where d— \z\ is the distance to the wall, while 

.1/2 



i?o = T^i'djv , = djRo , (6.4) 



i.e., Rq is the wall-stress Reynolds number and is the wall-length unit. By eliminating 
Rq in (6.3) we can write 

with some function h of two variables. Assuming that /i(0, oo) exists and noticing that 

, d—\z\\ , / (i— 

we obtain (as long as \z\ < do) that, for Rq large enough, the ratio 

- = -^/.(0,oo) 

is independent of Rq. This heuristic prediction will be confirmed later in a more rigorous 
way. 
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Finally, let us note that due to the symmetry of the physical setting, we can also 
assume that 

0-3(3^, y, -z, t; oj) = -as{x, y, z, t; uj) 

and therefore 

P{-z, t) = -/3{z, t) , a{-z, t) = a{z, t) . (6.6) 

7. Realizability conditions . Recall that the statistics of cr are subjected to the condition 
(2.15). In the present case this takes the form 

= {{didsai + d2dsa2)as) + {{diai + 82^2) ■ d^a^) + 2^3«^ - ^a/^ 
~ (9i<j2 • 92(72 - d2ai ■ 8102) + ((9i<Ji + 920-2) • Sacrs) - ((93(12 ■ d2a^) + {dia^, ■ dsai)) , 

whence 

^{ay -(3' ^ {diai ■ 92(72 - 92(71 ■ 9l(72) . (7.1) 

The meaning of a forces 

—d — z < (Ts{x, y, z,t;u>) < d — z for |^| < cZ . (7-2) 
In this case one can prove that the following conditions hold 

-d-z< (3{z) <d-z , a{zf < d'^ - z'^ - 2z(3{z) for \z\ < d . (7.3) 
Indeed, if P = P^^t denotes the probability distribution of as{z,t;u!) and 

p+=[ asP{das), p-=[ \as\P{das), 

J{a3>0} J{a3<0} 

then 

/? = {as) =/?+-/?-, r <{d + z)P{{as < 0}) , p+ < {d - z)P{{a3 > 0}) . 
Thus, 

{d + z)-^(3- + {d-z)-^(3+ < 1 , 

so that 

2dp_ <d^-z^- id + z)P . 
10 



On the other hand, 

= (^i) = / (^iPidas) <{d + z)p- + {d- z)P+ < 

< 2dp- + id- z)/3 <d^ -z^ - 2z/3 . 

This estabhshes the second inequality in (7.3). The first one is obvious. 

The Cauchy-Schwarz inequafity produces the supplementary constraint 

\P{z)\ < a{z) iov \z\ < d . (7.4) 

It is easy to check that the conditions (7.2) and (7.4) are also sufficient for the existence of a 
random variable (73{x, y, z; u>) satisfying (6.6) and (7.3) and statistically depending only on 
z. For any such (73, choose some homogeneous random vector [ai{x,y),a2{x,y)] such that 
^= {d^a^-d2a2-d2(T^-dia2)^0. Set = {2^)-^{{ay -2p')al Then a = (ai, (73, (73) 
has all the required statistical properties. We conclude that the inequalities (7.3) and (7.4) 
are the realizability conditions for the lengths a and /? in the VCHE (6.1). 

8. Comparing VCHE with the Revnolds equation . Comparing (5.1) and (6.1), we identify 
counterparts as, 

U^U , d,{wu) = u[{a^U'y" - iP'U)"] + po , 
d,{wv) = 0, V(P+(«;2)) = V(7r-poa;), (8.1) 

for a constant po. This identification gives 

{wv){z) = , 

-{wu)iz) = -poz - iy[{a^U'y'{z) - W'Uy{z)] , (8.2) 

and leaves (w^) undetermined up to an arbitrary function of z. A closure relation for 
— {wu) involving the third derivative U"'{z) also appears in Yoshizawa,^^ cf. equation (8) 
of Wei and Willmarth.^ 

^From (6.1) it follows that dxir = 7r2 is constant. Therefore integrating twice in z, the 
first equation in (6.1) gives 

- (3'{z))U{z) + v{a''{z)U'{z)y = ttq + ttiz - ^TTa^^ (8.3) 

with constants 7ri(z = 0, 1, 2). But the left hand side of (8.3) is symmetric under the change 
z I— > —2;, so TTi = and we obtain the following relation among the profiles of (5{z), a{z) 
and U{z), 

- P'{z))U{z) + u{a{zfU'{z)y = TTo - ^tts^^ for \z\ < d . (8.4) 
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For \z\ < do,P{z) = 0, a{z) = ckq > and (8.4) becomes 



-U{z) + alU"{z) = -ttq - ^TT2Z^ for \z\ < do . 



(8.5) 



Since U is symmetric in we obtain 



U{z) = a(l J ' % +6 1- 

\ cosh(do/Q!o)/ V % 



+ c for \z\ < do , 



(8.6) 



where the constants a, b and c satisfy the conditions 



c = U (±(io) , ttqI' = —a — 6(1 + 2q;o/c?o) ~ ^ , 772^' = —2h/dQ 



(8.7) 



It is worth mentioning here that with an antisymmetry condition for U{z) and with (8.6) 
changed accordingly, one may address turbulent shear flows (Couette flows) by the same 
analysis as developed in this paper. 
Integrating (8.4) on [— d, 0] gives 



.0 -, 

-V / {U{z) + P{z)U'{z))dz - a{-dfTo = Trod - -tts^^ 
J-d 6 

where we used (5.2) as well as £^'(0) = 0,/3(0) = and U{—d) = 0. Denoting 



(8.8) 



/•—do 

d 



U{z)dz + 



all ^ tanh — 



do 



ao 



-6 + c 



do_ 
d 



(8.9) 



allows (8.8) to be written also as 

-d 



/do ^ 
j3{z)U'{z)dz - ai-df'TQ = ttqcZ 7r2(i^ . 
-d 6 



(8.10) 



9. Empirical qualitative properties . It is universally accepted that the maximum of U is 
at 2 = (i.e. the center of the channel) and that U' {z) ■ z < Q ioi Q < \z\ < d. Also all 
experimental data show that U"{z) < over most of the channel. Thus 



i? := -t/av = 7^ / U{z)dz <Rc= -U{0) 



(9.1) 



12 



and (using the concavity property of U) 



R>- —^U{Q)dz=-Rc. (9.2) 

I' .1 —rl d Z 



Then (9.1), (9.2) can be given the form 



2 Rq tt* 

All the empirical evidence shows that 



1 1/(0) ^ B < CM . (9.3) 



R R 



2 < K — for i?o > 1 . (9.4) 



Throughout, the properties (9.3) and (9.4) will be taken as granted. 

10. The wall units representation . In the lower half of the channel, the mean velocity U 
can be expressed in wall units using the notation (l){rf) = U{z)/u^,r] = {z + d)/£^, with 
£^ = ly/u^ — d/Ro- In this representation, (8.6) becomes 

cosh^(l — T]/ Rq) 



a 






It* 








+ -I 









\ cosh^(l - ?7o/-Ro) 
1 - r//i?o 



1 - ?7o/-Ro 



+ </'(?7o), (10.1) 



for rjQ < rj < Rq, where ^ = d/a and r]Q = {d — c/q)/^* ~ cto/£* = Ro/C 
The definition of (p, implies R — J^" (l){r])dr]. Hence (10.1) gives 



aiRo-Vo) tanhe(l-go) ^ , 2b(i?o - r/p) , ^ ^ Ta., 

R= 1 ^ — + — i + <P{m){Ro - m) + (pi^jdr] . 

u* V ^(1 - Qo) J Su^ Jq 

To conclude this computation it is sufficient to approximate on (0, ?7o) by the piecewise 
linear function equal to ?7 for < < and (po + (r/ — t]o)<P'o for rj^ < rj < rjo, where 
4>o = Hvo), 4>'o = <P'{Vo) and ry* = {cpo - ?7o(^'o)/(l - </>'o)- We obtain 

R^ _ a(l - go) A _ tanhe(l-Qo) \ ^ 26(1 - gp) 
Rq V C(l - Qo) J 

+(1 - go)0o + (1 - 0'o)-^ (00^0 - iMR±^^ (10.2) 

where 

0'o = (a/«*)(C/^o) tanh(e(l - qo)) + 2(fe/«,)/i?o(l - qo) . (10.3) 



13 



Using this and solving for 0o gives an explicit function 0o = 0o(9o; R-, Ro'i o/m*, 
b/u^;^), namely 

Mqo; R, Ro; a/u,, b/u,;^ = Roi^B- ^ - i-[2(i?/i?o - C) + 9gi?o0(,]| (10.4) 

where 0'o is given by (10.3), 

5=(l-9o)(l-0'o) + 9o 

^ ^ a(l - go) r _ tanhe(l-Qo) \ ^ 2 6(1 - gp) 
tt* V C(l - go) / 3 

and the choice of the root (po in (10.2), (10.3) will be justified at the end of Section 11. 
Thus (10.1) becomes 

_co*|(l-^A (10.6) 

V cosh^(l - go) J 

b ( fl — ri/Ro\\ a b 



+ — 1 - " + 0o(go; i?, Rq\ — , — ; e) for goi^o < V < Ro ■ 

y \ 1 - qo J J 

In (10.6) the constants a/u^,,b/u^,^ and go may depend on Rq. As we will show below. 
Nature seems to choose them as constants (at least for large Rq). Recall that in Section 6 
we already gave a heuristic argument that ^ — d/a should be independent of Rq if -Ro (or 
R) is large enough. 



11. The off wall region . The empirical data up to now suggest that for a fixed channel 
there is a range [zi, Z2) (with 21^2 > 0) inside the channel such that for z in that range, 
the von Karman log-law is a good approximation to U{z), at least for R (or Rq) large 
enough. Since for those z we have 

C7(.2)-C/(.) = -ln^ = l(ln^-ln^) 
K z K \ a a/ 

(where k ~ .4 is the von Karman constant), U{z2) — U{z) is a function of z/d only (i.e. 
independent of Rq). We will posit now the following weaker condition. 

For R (or Rq) large enough, there exists a fixed range {zi/d,Z2/d) such that for z/d 
in that range, U{z2) — U{z) is a function of z/d, independent of Rq. 

Note that we make no assumption on the length of the range. The classical "defect 
law" of Izakson, Millikan and von Mises^-*^ (pp. 186-188) is the particular case of our 
condition when one of Zi^s is 0, and the range is assumed to be wide. 
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Passing to the wall units representation we can formulate our assumption as: There 
exists < qi < q2 < I, such that for QiRq < r] < q2Ro, (f){r]2) — (j){ri) is a function of 
q = tj/Rq only. Since we expect qo in (10.6) to be quite small, we will take qo < qi- 

We will prove now that under the above conditions, there exist absolute constants 
a* , 6* and ^* such that 

a ~ a*tt* cosh ^*(1 - go) 5 6~6*tt*(l-go) and ^ = (i/(Q;o^*) (H-l) 

where a, 6, ^ and qo are as in (10.6). 

Indeed let / be the function defined by 

/(g) = 4>{Q2Ro) - 4>{qRo) for gi < g < gs • (11.2) 

Then since go < gi we have from (10.6) 

/(g) = ao[coshe(l - g) - coshe(l - gs)] + 6o[(l - - (1 - 92)'] (11.3) 

where 

ao = (a/'u*)/cosh^(l - go) , bo = {h/u^)/{l - qof . (11.4) 

Writing (11.3) for g = gi, we obtain 

^ /(gi) - ao[cosh^(l - gi) - cosh^(l - ga)] 
°" (l-gi)2-(l-g2)2 

Then (11.3) becomes 

aogi^, q) = h{q) for gi < g < g2 (11.5) 
where, with cq an absolute constant, 

g{C, g) = coshe(l - g) - coshe(l - gs) - ci[(l - qf - (1 - gs)'] 
/^(O = fiq) - co[(l - g)' - (1 - q2f] , (11.6) 

and ao, ci^ are parameters, constant in g but which may depend continuously on Rq. Note 
that 

9{^,qi) = Kqi) = (z = 1,2) , ^(^g) < for gi < g < gs , C > . (11.7) 
Thus (with g = (gi + g2)/2) 

ao = h{q)/(iq) , (11.8) 
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and 



9{^,q)Kq) = Kq)^^^Q) for qi<q<q2 



(11.9) 



If ^ = ^(-Ro) were not constant, then (11.9) would hold for ^ in an interval [^1,^2] with 
< ^1 < ^2- Differentiating (11.9) with respect to ^ gives 

for ^1 < ^ < ^2,qi < 9 < 92- Introducing ^ = ^(1 — qf), it follows that 

sinhC = /io(0 + ^i(OC' + /i2(OcoshC for < C < 6, ^(1 - ^2) < C < ^(1 - ^i) , 

where ho, hi, /i2 are explicit functions of ^ only. Clearly this is impossible. 

We conclude from this contradiction that there are absolute constants qo,a*,b^ and 
^* such that 



cosh^*(l — go) — cosh^* ( 1 



Ro 



(l-^o)^-(l-f 



+ 00 (90; R, Ro; a*cosh^*(l - go),6*(l - go)^;C*) , 



(11.10) 



for qoRo <r]<Ro, where the function 0o (see (10.4)) actually depends only on qo, Rq and 
R. 

The formula (11.10) can be also written as 



4>{qRo) = 01 (go; q) + 0o(go; R, Ro) for go < g < i , 



(11.11) 



where 



h{qo;q) = a*cosh^*(l - go) ( 1 - 



cosh^(l — q) 
cosh^(l - go) 



+ 6*(l-go)^ 1- 



1-g 

1 - go 



for go < g < 1 



0o(go; R,Ro) = 0o(go; R, -Ro;a*cosh^*(i - go),fc*(i - go)^;C*) 



(11.12) 



For i?o ^ 00 from (10.4) and (10.5) we have 



R 1 

0(go-Ro) = 0o(go; R,Ro) — -C+ -goCo , 



(11.13) 
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where C is defined in (10.5) and Cq = Rq4>'q is constant according to (10.3). 

We can now explain the choice of 0o in (10.4). The other possible choice was 

Ro |s + y52-i-[2(r/ro-C) + ?2^o] } 
which would have given for Rq large enough 

(j){qoRQ) ~ 2Rq 

and consequently (j){Ro) > Rq which is contrary to the established facts (9.3), (9.4). 

12. The mean velocitv profile in the channel . Comparing the profile given by formula 
(11.10) with an experimental mean velocity profile, enables us to obtain the values a*, 6* 
and ^* as well as the smallest acceptable value for qq. In Figure 1, we compare our 
formula with experimental data^ for the Reynolds numbers Rq equal to 714, 989, and 
1608. As these Reynolds numbers are small, a* and 6* have not reached their asymptotic 
values. It appears, however, that ^* has reached its asymptotic value. We therefore allow 
a* and 6* to vary slightly with Rq while holding ^* constant to fit the data. It turns 
out that ^* = 35 and = 1/^*. Note that this choice of q^ corresponds exactly to the 
condition that |ci — do] = a. 

13. The Revnolds shear stress . The shear Reynolds stress is —{uw) (see Section 5). Since 
(tttu) = 0, one must have 

-{uw){z) = -To^ - uU'(z) for \z\<d. (13.1) 

On the other hand U = U and —{uw) is also given by (8.2) with an appropriate constant 
Pq. For \z\ < do, since a{z) = ao = l/^*,P{z) = 0, (8.2) reduces to 

-{uw){z) = pQZ - valU"'{z) (13.2) 

and 

ValU"'{z) = VU'{Z) - TT2Z . 

Introducing this in (13.2) we see that (13.1) and (13.2) are compatible if pq is given by 

Po = --kq - To/d . (13.3) 
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Taking the wall units representation in (13.1) we obtain our theoretical Reynolds shear 
stress 



>"^=l-^-0'(r;) = 



To -Rq 

= 1 - ^ + ^^*sinh^, fi-A\+^fi-^\ for g,i2o < 77 < i?o . (13.4) 
no Kq \ no / Kq \ no y 

Figure 2 compares the corresponding experimental and theoretical Reynolds shear 
stresses. We use the same values for a*, 6* and ^* as before. The agreement in shear 
stresses does not extend as close to the wall as the mean velocity profiles did. The empirical 
matching of the mean velocity profiles as well as the Reynolds shear stresses are both given 
with = -\/3/^*. We note that the consistency of this closure and the experiments found 
in the trends followed by the Reynolds-stress profiles in Figure 2 is an exacting test of 
the fidelity of the mean velocity profiles as well as a test of the Reynolds stress relation 
predicted by equation (13.4). 

14. The near wall region . As already mentioned above, in the near wall regions (i.e. where 
< ?7 < r/o = Q*Ro and 2i?o — q*Ro < < 2i?o), P may be non-zero and a may depend 
(as does /?) on r] and Ro- The VCHE (8.4) in the wall units representation takes the form 

(1 - RoP'ivmv) - RU&ivf<p'{v)y = /o + 3/1 (^1 - 

for 0<r]<Ro (14.1) 
i.e., in the whole lower half of the channel. In (14.1) we used the notations 

a{ri) = a{z)/d , P{rj) = (5{z)/d (14.2) 
where d + z = ry^*,^* = d/Ro, and 



fo = — ttq/z^w* = a*cosh^*(l — g*) -|- 6* 



(1 - q.f + ^ 



-I- 0(g*-Ro) 



fi = T^id^jvu^ = -6*/3 (14.3) 
Of course we have 

oiiri) = , Piv) = for q^Ro < V < Ro , (14.4) 
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but the VCHE (14.1) does not define a{r]),P{r]) near the wall (i.e. for < 77 < q:^Ro). 
However, (14.1) gives some qualitative information on the behavior of a and $ in the near 
wall region. Indeed integrating (14.1) we obtain 



1 rRo rRo 

— / cf>in')dv' + / pw)cf>'W)dv' + mm + Rom^^' 



iv) 



for all < ?7 < Rq. Thus (using \P\ < 2, see(7.3)) we find 



which in turn goes to zero when Rq — 00, by virtue of (9.4). Thus, for rj fixed. 



a{'r]) (p'irj) — > for Rq ^ 00 . 
In particular, (using 0'(O) = 1 and (9.3)) we obtain 



5(0) ^ , p{0) ^0 for Rq^oo 



and, due to (14.4), 



0'(g*i?o) ^0 for i?o ^ 00 . 
Moreover, writing (14.5) for = 0, 

R ~ 

-W- + / mW{v')dr]' + Rom^ ^fo + fi 
-n-0 Jo 

and subtracting (14.5) from (14.9) we obtain 
1 rv 



R 



^0 



(14.5) 



(14.6) 



(14.7) 



(14.8) 



(14.9) 



m)dr/ + / m)<p'W)dr/ - mm + roHoy - m'<p\v)] = 



— /o ^ + /i 



Rq 



(14.10) 



Fixing r] and letting Rq — > 00, from (14.10) we infer that 

r m)<t>'{v')dv' - mm + roHoy - m^4>'{v)] ^ o . (m.h) 
^0 



19 



But for rj fixed such that (f)'{r]) remains away from when i?o we have P{rj) ^ 

because of (7.4). By virtue of Lebesgue's dominated convergence theorem we conclude 
that 

Ro[a(Of - a{r}f^'{r})] when i?o ^ oo (14.12) 

provided that (/>'(?7) stays away from 0. 

It is instructive to connect (14.12) with the functions g and h considered in Section 
6. First we recaU that 

g{Ro, v) = m , f0iv).v) = Hv) ■ (14.13) 
Assuming that P and P' are continuous across rjo = q*Ro allows us to conclude that 

g{Ro, Roq) = - J {q' - Q*)[-q^) 9{Ro, Roq')dq' for < q < q* . 

Since qf* is small, assuming that 

{g^)g{Ro,Roq') 

exists and is not zero, we obtain that 

9{Ro,v) - ll{Ro){q* - ri/Rof • (14.14) 

But /3(0) > 0, so 7(-Ro) > 0. If 7(-Ro) = 0, we can proceed in a similar way by involving a 
higher derivative of g at ryo from the left. In all cases we have ended with a representation 

g{Ro, 7]) ~ goiRo)gi{v/Ro) for < < q.Ro (14.15) 

where ^ 0, > and gi{q) is a decreasing function of q with gi{q^.) = g[{q*) = 0. 
From (14.15), (14.13) and (14.12) we now obtain 

h\go{Ro)gm,0) - h\go{Ro)gi{v/Ro),v)<l>'{v) ~ (14.16) 

for Rq large enough. In (14.16), gi{r]/RQ) — > 5^1(0) for Rq ^ 00 and rj fixed. These 
arguments suggest that 

h{p,rj)r^h{P,0)/,/^ (14.17) 

provided P is small and 1/0' (ry) is bounded when i?o — > 00. It is not clear if assuming 
equality in (14.17) is a judicious approximation of the function h. 
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A major difficulty in fine tuning our approach near the wall resides in the unavailability 
of experimental data in the near wall region for large Reynolds numbers. To test whether 
the VCHE (14.5) is still valid in this region we extrapolated the experimental profiles 
in Figure 1 into the near wall region according to Panton^ to obtain a from (13.5). For 
simplicity of the graph, we will display the a profile only for i?o = 1608, which is the highest 
Reynolds number in Figure 1. As illustrated in Figure 3, we find that the realizability 
conditions (in Section 7) are satisfied for appropriate choice of 7(-Ro) in (14.14) and (14.13). 
Clearly a lies between the upper constraint (7.3) and the lower constraint (7.4) in the near 
wall region. Thus, our basic ansatz is consistent with Panton's theory. 



15. Pipe flows and prediction . All the preceding considerations on turbulent channel flows 
can be suitably applied to turbulent pipe flows. The substantial difference between the 
mathematical treatment of the two types of flows, is that for pipes, the cosh function is 
replaced by the flrst modified Bessel function^^ 

w = E(;^(^) . (15.1) 

n=l ^ ' ^ ^ 

For instance the basic formula (10.11) becomes 
0(^) = — 

+ 4>{7]q) for f]o<V<Ro- (15.2) 

For pipe flows, experimental data for quite large Reynolds numbers are available (see 
Zagarola^). For these Reynolds numbers it is reasonable to assume that a*, K, and ^* 
have each reached their asymptotic values. In Figure 4 we compare our proflles with 
experimental data of Zagarola.^ We obtain the a*, 6*, ^* and by using the experimental 
data for R = 98, 812 and use the von Karman drag law, R/Rq ~ logi?o, to obtain proflles 
for R = 3, 089, 100 and 35, 259, 000. See also Chen et al.^^ for additional discussion and 
numerical details for these comparisons. 

We note that our predictions are consistent with the von Karman log law,^^ the 
Barenblatt-Chorin power law,^° as well as with the presence of the 'chevron' near the 
center of the flow. Our approach shows a logarithmic proflle for 0.02 Rq < r] < 0.2i?o 
and a chevron near the center of the channel. The Barenblatt-Chorin power law-^° may 
represent the transition in the proflle from the log law to the chevron. Although our 
approach is in good agreement with the experiments,^ we note that it has been argued 



lom-Vo/Ro) 



^ b f f l-v/Ro 
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that the experimental mean velocity profiles are too low for high Reynolds numbers. ° 
Finally, we observe that the chevron may reflect the fact that, on the attr actor of the 
dynamical system in the phase space of the turbulent flow, the Poiseuille-Hagen flow is 
recurrent. 
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Figure Captions 



FIG. 1. The mean velocity profile in the channel for the constant-a viscous Camassa-Holm 
equation compared with the experimental data of Wei and Willmarth.^ 

FIG. 2. Reynolds shear stress in the channel compared with the experimental data of Wei 
and Willmarth.^ 

FIG. 3. The statistical compatibility of the VCHE with the theory of Panton^ in the near 
wall region. 

FIG. 4. The mean velocity profile in the pipe for the constant-a viscous Gamassa-Holm 
equation compared with the experimental data of Zagarola.^ 
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